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Abstract 

The so-called renormalization group (RG) method is applied to derive kinetic and 
transport equations from the respective microscopic equations. The derived equa- 
tions include Boltzmann equation in classical mechanics, Fokker-Planck equation, a 
rate equation in a quantum field theoretical model. Utilizing the formulation of the 
RG method which elucidates the important role played by the choice of the initial 
conditions, the general structure and the underlying assumptions in the derivation 
Q \ of kinetic equations in the RG method is clarified. It is shown that the present 

formulation naturally leads to the choice for the initial value of the microscopic dis- 
tribution function at arbitrary time to to be on the averaged distribution function 
to be determined. The averaged distribution function may be thought as an integral 
constant of the solution of microscopic evolution equation; the RG equation gives 
the slow dynamics of the would-be initial constant, which is actually the kinetic 
equation governing the averaged distribution function. It is further shown that the 
averaging as given above gives rise to a coarse-graining of the time-derivative which 
is expressed with the initial time to, thereby leads to time- irreversible equations even 
from a time-reversible equation. It is shown that a further reduction of Boltzmann 
equation to fluid dynamical equations and the adiabatic elimination of fast vari- 
ables in Fokker-Planck equation are also performed in a unified way in the present 
method. 
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1 Introduction 



X 

Statistical physics, especially of non-equilibrium phenomena may be said to be a col- 
lection of theories on how to reduce the dynamics of many-body systems to ones with 
fewer variables]!]]. BBGKY( Bogoliubov-Born-Green-Kirkwood-Yvon) hierarchy which 
is equivalent to Liouville equation hence time-reversible can be reduced to the time- 
irreversible Boltzmann equation H given solely in terms of the single-particle distribution 
function for dilute gas systems^. The derivation of Boltzmann equation by BogoliubovQ] 
shows that the dilute-gas dynamics as a dynamical system with many-degrees of freedom 
has an attractive manifold |J spanned by the one-particle distribution function, which is 
also an invariant manifold |j| . Boltzmann equation in turn can be further reduced to the 
hydrodynamic equation (Navier-Stokes equation) by a perturbation theory like Chapman- 
Enskog method^] or Bogoliubov's method[[|, 0. Langevin equation which may be time- 
irreversible can be reduced to the time-irreversible Fokker-Planck equation with a longer 
time scale than the scale in Langevin equation||. Two basic ingredients are commonly 
seen in the reduction of dynamics, which are interrelated with but relatively indepen- 
dent of each other: (i) The reduced dynamics is characterized with a longer time scale 
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than that appearing in the original (microscopic) evolution equation, (ii) The reduced 
dynamics is described by a time-irreversible equation even when the original microscopic 
equation is time-reversible. The fundamental problem in the theory of deriving kinetic or 
transport equations is to clarify the mechanism of and to implement the above two basic 
ingredients. 

It seems that the basic notions to implement the two basic ingredients (i) and (ii) are 
given by (1) the coarse-grained time- derivative and (2) the choice of the initial conditions 
in solving the microscopic equations, respectively: 

(1) In an attempt to characterize hydro-dynamical processes microscopically, Mori pointed 
out that time derivatives appearing in equations which define transport coefficients are 
"the average" of time derivatives describing microscopic dynamics |§. His definition of 
the macroscopic derivative of an observable F is 



A(F)(t) = -{(F)(t + r) - (F)(t)} = - [ T ds^-(F)(t + s), (1.1) 
01 T r Jo as 



where r is some time scale between microscopic (mean free) time and macroscopic (re- 
laxation) time. An important point is that r is finite. The idea of the coarse-graining of 



time in kinetic and transport equations were first given by Kirkwood||10|; see also ||TT| for 
a rigorous formulation. 

(2) The importance of the choice of the initial condition in the derivation of kinetic equa- 



tion is noticed and emphasized in the literatureS H, 02, 13, 14]. For instance, Kawasaki 



clarifies in an excellent monograph ||14|| that the initial value of the microscopic distribu- 



tion function before averaging must be given solely in terms of the averaged distribution 
to obtain a closed equation for the distribution function of macroscopic slow variables, 
which is equivalent to the construction of an invariant manifold mentioned above. He also 
clarifies that by this initial condition, the dominating class of states ( "typical states" ) are 
selected which leads to an increase of the entropy, while exceptional states from which the 
entropy would increase could be unstabilized to the dominating typical states by a mech- 
anism producing chaotic behavior. Bogoliubov [|J also gave exactly the same scenario as 
described above in his asymptotic theory for the reduction to Boltzmann equation; see 
chap. 9 of [§. 

In this paper, we apply the so-called renormalization group (RG) ]15| method to derive 
and reduce kinetic equations to a slower dynamics JIB], [T3|, [B], |2"0fl : We show that if the 



RG method is properly formulated so that the essential role played by the choice of initial 
conditions is manifest as is done in Jl9|, |20|, the general structure in the derivation of 
kinetic equations in the RG method is exactly in accordance with the basic principles (1) 
and (2) mentioned above. 



The reduction-theoretical aspect of the RG method [|T^] and the improved formulation 



given in |l9j were reformulated mathematically with the notion of invariant manifolds fa- 



miliar in the theory of dynamical systems || by Ei, Fujii and one of the present author 



In [20], it was shown that the perturbative RG method can be used to construct invariant 



manifolds successively as the initial values of evolution equations using the Wilsonian 



RG |T5|, [21], p2| ; the would-be integral constants, which have one-to-one correspondence 
with the initial values, in the unperturbed solution, constitute natural coordinates of the 
invariant manifold. It was also shown that the RG equation determines the slow motion 
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of the would-be integral constants on the invariant manifold of the dynamical system, 
hence a reduction of evolution equation is achieved. 

We shall show that the straightforward application of the RG method as formulated 
in fl"9| , [2C|j naturally leads to the choice for the initial value of the microscopic distribution 



function at an arbitrary time to to be on the averaged distribution, which is an implemen- 
tation of (2) in the RG method, thereby leads to time-irreversible equations even from a 
time-reversible equation. The averaged distribution function may be thought as an inte- 
gral constant of the solution of microscopic evolution equation. The RG equation gives 
the slow dynamics of the would-be initial constant, which is actually the kinetic equation 
governing the averaged distribution function. It will be further shown that the averaging 
as given above automatically gives rise to a coarse-graining of the time-derivative, which 
is expressed with the initial time to. This shows that the initial time to has a macro- 
scopic nature in contrast to the time t appearing in the microscopic equations, which is 
an implementation of (1) in our method. 

It should be noticed here that the RG equation has been already applied to kinetic 
equations both in classical and quantum regimes by others^, ^4|: In p3"] , it was as- 



certained that Boltzmann equation is a renormalization group equation on the basis of 
the work by M.S. Green for the uniform system which shows that the perturbative 
solution for BBGKY hierarchy exhibit a secular term, and a sketch was given to derive 
Fokker-Planck equation from a simple Langevin equation noticing again an appearance 
of a secular term. On the basis of these two examples, they claimed that all other kinetic 
equations are also RG equations. Boyanovsky, de Vega and their collaborators ]24| have 
derived a Boltzmann equations for quantum field theories including gauge theories in the 
prescription of Illinois group [1£|. They noticed that the existence of a mesoscopic scale 
which is microscopically large but macroscopically small is essential for the applicability 
of the RG method to derive kinetic equations as well as the validity of the description in 
terms of a kinetic equation itself; we remark here that the existence of such mesoscopic 
scales may be considered as the so-called intermediate asymptotics f26fl , which was actu- 



ally the basic observation for the application of the RG equation to asymptotic analysis of 
differential equations [15]. A model dealt in |23|, |24| will be retreated in our formulation, 
and implicit assumptions in their treatment will be made explicit so that the roles of the 
initial conditions and the scale transformation of the time-derivative will become clear for 
the RG method to lead to kinetic or transport equations. 

In section 2, we shall deal with Langevin equation and derive the Fokker-Planck equa- 
tion as a typical problem of dynamical reduction leading to a kinetic equation in the RG 
method. We shall summarize the basic structure of the reduction given by the RG method. 
One will see that a similar definition to Eq.(l.l) of the macroscopic time-derivative nat- 
urally emerges in the RG method. In section 3, Boltzmann equation is derived from 
Liouville equation of the classical mechanics; we shall clarify the difference between the 
present method and the one by Bogoliubov. In section 4, we apply the RG method to 
derive a Boltzmann equation(rate equation) in quantum field theory in a way where it 
is transparent how the initial value are chosen and the coarse-grained time-derivative is 
introduced in the RG method. This section is a recapitulation of the work by Boyanovsky 
et al in our point of view. In section 5, the RG method as formulated in |20| is applied 
to obtain the fluid dynamical limit of Boltzmann equation. This is an example to reduce 



3 



a kinetic equation to a further slower dynamics, which appear quite often in statistical 
physics reflecting the hierarchy of the space-time of the nature. Another typical problem 
in this category is to derive Smoluchowski equation from Kramers equation. In section 6, 
we show that it is possible to develop a systematic theory based on the RG method for the 
adiabatic elimination of fast variables in Fokker-Planck equations. Section 7 is devoted to 
a summary and concluding remarks. In Appendix A, the calculational procedure sketched 
in [23] for deriving Fokker-Planck equation in the RG method is properly elaborated and 
worked out so that the essential role of the initial condition in the procedure is fully rec- 
ognized; it is shown that the initial value to obtain the "stochastic distribution function" 
must be actually the averaged one for the RG equation to give Fokker-Planck equation in 
conformity with the observation by Kawasaki and others mentioned above. 
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2 From Langevin to Fokker-Planck equation 



In this section, we apply the RG equation to derive the Fokker-Planck(FP) equation|| 
starting from the stochastic Liouville equation [27|| corresponding to Langevin equation^. 
This derivation is thought to be a typical one for the reduction of evolution equations 
appearing in non-equilibrium physics ||14|| . It will be clarified that the initial values of the 
stochastic distribution function at arbitrary time to are naturally chosen to be on the 
averaged distribution function for the RG equation to derive the FP equation governing 
the averaged distribution function. We shall also notice that the time derivative in the 
RG equation which will be converted to the derivative in the FP equation is with respect 
to a macroscopic time, hence the coarse-graining of time is automatically built in in the 
present RG method. 



2.1 Application of the RG method to a generic equation Langevin 
equation 

We consider the following generic Langevin equation with Ri (i = 1, 2, n) being stochas- 
tic variables; 

fin 

^ = h{u) + g(u)R, (2.1) 

where u = t (u\, u 2 , u n ), h = t (hi, h 2 , h n ), g a n times n matrix and {g{u)R) i = 
J2j9ijRj- Notice that the noise enters multiplicatively. Here we assume without loss of 
generality that the noise has the vanishing average, 

(R(t)) = 0, (2.2) 

where (0(t)) denotes the average of 0(t) with respect to the noise R. Let f(u,t) be the 
distribution function with R(t) given; the continuity equation reads 

^ + V u -W(«,t)) = 0, (2.3) 

where V u = J2id/dui and v = du/dt is the velocity of u, which is given in ( |2.1| ). 
Inserting ( |2.1|) into (|2.3| ), one has the Kubo's stochastic Liouville equation [f27f , 



dt 



-V u -[(h + gR)f). (2.4) 



Whereas it is rather easy task to derive the FP equation in an exact way if the noise is 
Gaussian, it is formidably difficult if the noise is non- Gaussian^. Although the present 
approach is based on the perturbation theory and of approximate nature, it will be found 
that the first order calculation suffices to derive the exact FP equation when the noise 
in Gaussian, and furthermore that the method is applicable even to non-Gaussian noises 
without difficulties. 

The solution to (|2.4j ) with the initial condition given at t = to is formally given by [p7f 
f(u,t;t ) =Texp[/ dsL(s)}f(u,t ;t ), (2.5) 

Jtn 
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where 

L(s) = -V u -(h(u)+gR(s)), (2.6) 

and T denotes the time ordering operator. The initial distribution f(u,t ;t ) will be 
specified later and found to play a significant role in the present method. 

Now we are interested in the averaged distribution function P(u, t; to) which is defined 
as an average of /(u, t;to) with respect to the noise R, i.e., 

POM; to) = (Texp[/'d S L( S )]/(u,t )>. (2.7) 

We take an interaction picture dividing the "Hamiltonian" L as follows; 

L = L + Li, (2.8) 
L = -V u h, L l = -V u gR. (2.9) 

We first define U (t) as the time-evolution operator governed by the unperturbed "Hamil- 
tonian" L , 



U (t,t ) = Texp[/ t dsL (s)], (2.10) 

Jto 



ft 

which satisfies the evolution equation 

J^o(Mo) = L U (t;t ), (2.11) 

with the initial condition 

U (t ,t ) = 1. (2.12) 

Then to incorporate Li, we define another microscopic distribution function pi(u, t;to) 
by 

f(u,t;t ) = U (t,t ) Pl (u,t;t ). (2.13) 

We remark that the initial values of / and p\ at t — t coincides with each other, which 
we take to be equal to the averaged distribution function P(u,t ) at t — to as depicted 
in Fig.l; 

f(u,t = t ;t ) = Pl (u,t = t ;t ) = P(u,t ). (2.14) 

One will recognize that this choice of the initial condition is inevitable for the RG equation 
to be identified with Fokker-Planck equation. 

One can easily verify that pi(u, t; to) is formally solved to be 

p 1 (u,t;t ) = Texp[[ dsL^s; t )]pi(w, t ; to), (2.15) 
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Distribution functions 




Figure 1: The initial value of f(u,t = t ;t ) = p±(u,t = to; to) at t — to is set on 
the averaged distribution function P(u,to), which will be found to obey Fokker-Planck 
equation as the RG equation. 



where 



C 1 (t;t )=Uo L (t,t )L 1 (t)Uo(t,to), 



is an "interaction Hamiltonian" in the "interaction picture" . 
Thus we obtains the compact form of P(u, t; t ) as follows, 

P(u,t;t ) = (U (t,t )p 1 (u,t;t )), 

= U (t,to)(Texp[f ds£i(s;* )])P(tMo), 

Jt 

= U (t,t )S(t;t )P(u,t ), 
where we have used the fact that pi(u,t = to, to) = P(u,t ) and the abbreviation 

S(t;t ) = (Texpif 1 dsd^to)]). 

J to 

The computation may be performed in a perturbative way: 

S(t; t ) = 1 + T f dsiC^s)) + \t f ds x f ds 2 (C 1 (s 1 )C 1 (s 2 )) + ... 
J Iq Z J to J to 

= 1 + \t f d Sl f ds 2 T( Sl ,s 2 ) + ... 



[2.16) 



(2.17) 
(2.18) 
(2.19) 



where we have put 



T(si,s 2 ) = (£i(si)£i(s 2 ))- 



(2.20) 
(2.21) 
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If the noise is stationary, which we shall assume from now, T(si, S2) will be a function of the 
difference si — S2', furthermore, owing to the time-reversible invariance of the microscopic 
law, T(si, s 2 ) becomes a function of the absolute value \si — s 2 \, i.e., T(si, s 2 ) = r(|si — s 2 |). 
Then one has for t > to, 

S(t; t ) = 1 + (t - t )G(t - t ) + • • • , (2.22) 
where we have put for t > 

G(t) = f' dsT(s). (2.23) 
Jo 

If we stop at the second order approximation, we have 

P{u, t; 1 ) = U{t- to) [l + {t- t )G(t - t )]P{u, t ). (2.24) 

Notice the appearance of the secular term which indicates that the above formula is only 
valid for t around to- 

Now we apply the RG equation to ( |2.24|) which reads 

dP(u,t;t )/dt \ to=t = 0, 

which leads to 

d to Uo(t,to)\ to=t P(u,t) + d t P(u,t) - G(0)P(u,t) = 0, 

where use has been made that f/o(t ,t ) = l,Vt . Noticing that d to Uo(t,to)\t =t = —Lo = 
V u ■ h, we arrive at Fokker-Planck equation, 

d t P(u, t) = -V u ■ hP{u, t) + G(0)P{u, t). (2.25) 

The concrete form of G(0) depends on the character of the noise R(t). This is one of the 
main results of this section. 

To see that (|2.25|) is the desired equation, let us evaluate G(0) for a simple Gaussian 
noise given by 

(RityRjit')) = 2<%A6(t - *')■ ( 2 - 26 ) 

For this case, one has 

r(s) = Uo l d i gi j d k gki^D j 5 j i5(s), 

where di = d/dui. Then G{t) is evaluated as follows; 

r f 1 
G(t) = ^ dsT(s) = -Uo l d i g ij d k g k i2D j 8j h 

= G(0). (2.27) 

Here we have used the identity 9(0) = 1/2, in accordance with the Stratonovich scheme 0. 
Notice that G(t) in this case is independent of t. 

Inserting G(0) thus obtained into fl2.25|) , one has the familiar form of Fokker-Planck 
equation for the multiplicative Gaussian noise, 

d t P(u, t) = -V u ■ hP(u, t) + DjdigijdkgkjPiu, t). (2.28) 

This shows that the initial distribution P(u, to) satisfies Fokker-Planck equation and 
justifies the identification of the initial distribution with the averaged one made in Eq. 
( P-14| ), as depicted in Fig. 1. 
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2.2 Discussion 



Firstly, it is noteworthy that we have been naturally led to identify the initial values of 
the microscopic distribution function f(u, to, to) before averaging with the averaged value 
P(u,to) at an arbitrary initial time t = to- As mentioned in §1, the necessity to take 
such an initial condition to achieve reduction of evolution equation was advocated by 
Bogoliubov[[|] and others[13, 14] including Boltzmann[Q. Secondly, this means that the 



nature of the initial time to in the RG method is completely different from that of the 
time t in the stochastic equation (microscopic equation); t represents the coarse-grained 
time describing the variation of the averaged quantity, and the derivative dt in the RG 
equation is a macroscopic time-derivative. Again as mentioned in §1, this coarse-graining 



of time was also noted by others |T0|, || |TT| in different approaches. 

This automatic averaging and the appearance of the macroscopic time-derivative given 
in the RG method may be generically understood as a generalization of the scheme given 
in §2 of | 20|| : First discretize the variable u — > Ui and write as P(u,t)(ui,t) = Xj(t) and 



use a vector notation X = (Xx, X 2 , ■■■)■ Thus the discretized stochastic Liouville equation 
with the initial value X(to) at an arbitrary time to may be solved perturbatively, and the 
solution is denoted as X(t; to, X(t )), which satisfies the initial condition 

X(t ;t ,X(t Q )) = X(t ). (2.29) 

We could solve the same equation with the initial condition given at a shifted initial time 
t = t + At; 

X(t + At; t + At, X{t + At)) = X{t + At). (2.30) 

We suppose that the time difference At is macroscopically small but microscopically so 
large that it may be taken as infinity. For the time t between to and to + At, i.e., 
to < t < to + At, the perturbation should be valid. If t — to and At — > oo in the 
microscopic scale, we may anticipate that the system is relaxed to the averaged trajectory 
X(t) and have 

X(t; to + At, X(t + At)) ~ X(t; t , X(t )), 
which implies that the macroscopic time derivative S/5to vanishes, 

6X _ X(t; tp + At, Xjtp + At)) - X(t; t , Xjtp)) 
5t At ' 1 ' } 

dX dX dX 
~ dto~ l °= t + dX'Ho~' ( J 

Notice that in the macroscopic scale, the equality t — t ~ to + At should be taken for 
granted. This is the RG equation underlying the derivation of Fokker-Planck equation 
and also other transport equations including kinetic equations as will be shown in the 
proceeding sections. 

In the following sections, we shall show that the general structure elucidated in this 
section persists in other important examples of the reduction to kinetic equations. 
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3 From BBGKY hierarchy to Boltzmann equation 



In this section, we apply the RG method to derive Boltzmann equation for a classical 
dilute gas starting from BBGKY hierarchy 0. As is well known, Bogoliubov first derived 
Boltzmann equation from BBGKY hierarchy in his classic paper ||. His derivation starts 
from an ansatz that the many particle distribution function depends on time only through 
the one-particle distribution function and uses a special perturbative expansion method. 
His approach is actually an application and generalization of the asymptotic theory by 
Krylov and Bogoliubov (KB) successful to non-linear oscillators 0. Here we do not make 
that ansatz and start from the naive perturbation theory. We will see how the ansatz given 
by KB can be incorporated in the RG method. The importance of the initial condition 



again emerges, which is in accordance with the observation given in [28[ . This implies that 
the appearance of a secular term|23| does not constitutes the final story for the derivation 
of Boltzmann equation. 

3.1 Derivation of Boltzmann equation in the RG method 

Consider a system of N identical classical particles enclosed in a volume V. Following 



the notation of p9|, we denote the 2-th particle's phase space coordinate by Xi = (r i ,p i 



The Hamiltonian of the system is 

ff = E^ + 5E^l-.-ol)- (3-D 

i=i i^j 

We assume that the potential U depends only on the relative distance of two particles and 
that its range d is much shorter than the mean free path I. The iV-particle distribution 
function /jv(xi, • • • , x$r, t) is normalized as 

f i +\ 0£i drid Pi 1 to o\ 

f N (x u --- J x Nj t) — = 1. (3.2) 

We define the s-particle distribution function by 

f s (x 1 , ■ ■ ■ x s , t) = J f N (x 1 , ■ ■ ■ ,x N , t) dXs ^ _ ^ N . (3.3) 

Then the normalization condition for f s becomes 

r AH 

J fs{xi, ■ ■ • , x a , t)dxi ■■■dx s = ^ N _ s ^ - N ^ ( 3 - 4 ) 

from which we see that f s is of s-th order in the particle density n = y. We assume that 
n < 1. 

The kinetic equation for f s is obtained by integrating Liouville equation 4/jv = over 
x s+ i, • • • , Xjv- Equations for fi and / 2 read 

-j t h{xi,t) = (^ + iL° 1 )f 1 (x 1 ,t) = - J dx 2 L' 12 f 2 (x 1 ,x 2 ,t), (3.5) 
d r 

dj.M x ii x 2, t) = (— + iLi 2 )f 2 (xi, x 2 , t) = - J dx 3 (iL' 13 + iL'^f^x^ x 2 , x 3 , t)(3.6) 
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where 



12- 



i? = -^'f. ^ =j «2^£i.( » _ » ). (3.7) 



These are the first two equations of the BBGKY hierarchy which is a series of equations 
relating the evolution of f s to A+i- Our goal is to derive an equation (or equations) 
which captures the essence of the system's dynamics described by BBGKY hierarchy. 
In the language of the the theory of dynamical systems [H, we wish to construct a low- 
dimensional invariant manifold in the (practically) infinite-dimensional functional space 
spanned by {f s } and derive the reduced equations of motion on it. 

Whereas Liouville equation or, equivalently, BBGKY hierarchy describes microscopic 
collisions between particles in detail, what interests us is the macroscopic variation of the 
system caused by the accumulation of many collisions. More concretely, we wish to know 
the variation of the system over the space-time scale much longer than the collision radius 
and the collision time and much shorter than the mean free path and the mean free time. 
Such scale is called the mesoscale. The derivatives appearing in ( |3.5| ) and (|3l6|) are, so to 
speak, microscopic derivatives, while those appearing in kinetic equations are macroscopic 
derivatives. We must take into account their difference when deriving kinetic equations. 

Following [p0|1 , suppose that we have found the solution to BBGKY hierarchy {/ s (, t)} 
up to an arbitrary time to- With the initial condition {AG to)} we t r Y to solve (|3.5|) and 
( p.6| ) by the perturbative expansion in the density (virial expansion) to obtain a solution 
f s (t;to) around t ~ to- Recalling that f s is of s-th order in the density, we expand as 
follows. 

fi(xi,t) = ff(x 1 ,t)+fl(x 1 ,t)+f*(x 1 ,t) + ---, (3.8) 

f 2 (xi,X2,t) = f°(xi, x 2 , t) + fl{x\, x 2: t) -\ , (3.9) 

f 3 (x u X2,x 3 ,t) = f^(xx,x 2 , x 3 , t) H , (3.10) 

where fl{x\, ■ ■ ■ , x iy t) is of (i+j)-th order in the density. Substituting the above expansion 



in (|3.5|) and (|3.6|) , we get 



f t f° l (x 1 ,t) = 0, (3.11) 

j t ft(x 1 ,x 2 ,t) = Q, (3.12) 

+ - • 7T-)/lOM) = / d ^^-U(\ ri -r 2 \). JLf2( Xl , X2 ,t), (3.13) 
at m or i J or\ op 1 

where we have dropped terms which result in the surface integral. We also expand the 
initial condition 

A (a*, to) = fiixt^o) + fl(x u t ) H , 

f 2 (x u x 2 ,t ) = f 2 (xi, x 2 , t ) H • (3.14) 
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Equation ( |3.11| ) and ( |3 . 1 2| ) are easily integrated: 



f%(x u x 2 ,t) 



iL 12 {t WfO^Xi, X 2 , to) = f 2 ( x 10i x 20, to) 



where 



%iO\ x li x 2, t, to) 



1,2 



(3.15) 



(3.16) 



are positions and momenta of the particles 1 and 2 at time t under the influence of the 
2-body Hamiltonian 



2m 2m 



H' 21 = -^ + -^- + Ul\r -r 2 \). (3.17) 
The initial values fi(xi,t ) and f 2 (xi,x 2 ,t ) may be considered as the integration con- 



stants of the lowest-order equation. In the RG method as formulated in JT9|, BO], the 



integration constants will constitute the coordinates of the zeroth invariant manifold | 
The decisive step of the present approach is to choose the initial condition as follows 



f 2 (x!,X 2 ,t ) = fi(x 1 ,to)fi(x 2 ,t ) 



(3.18) 



irrespective of the distance between r% and r 2 . The underlying picture of this choice is 
that the system is so dilute that the two particles at an arbitrary time to are most probably 
located at distance much longer than the collision radius d, so that the correlation of the 
two particles is negligible and f 2 can be set to the product of one-particle distribution 
functions. We remark that a probabilistic nature enters at this point [|30|j . 

The integration of (|3.13| ) from to to t with - t — 1 (v is the average velocity), which 
implies that t — 1 is small in the macroscopic scale, gives 



fl(x u t) 



e - lL 0(t-t ) f l {xuto) 

+ fdt'e-^^ [ dx 2 ^-U{\ ri - r 2 \) ■ ^/ 1 (x' 10 ,t )/ 1 (x / 20) toX3-19) 











dp{ 



where we have used (|3.15|) and (|3.18|) , and x' 10 and x 20 are given by (glBD with the 
replacement t — > t' . We remark that the condition (- ^> t — t ) is also required for the 
expansion in the density to be valid J3J. In (|3.19|) , only r 2 for \r% — r 2 \ < d contributes 
to the integral. In this region, we can write 



r i0 ~ r i 



Pa 

m 



it' -to). 



(3.20) 



for a microscopically large period if — to 3> -. Here we have neglected vectors whose 
magnitudes are of order d. Then the perturbative solution in the mesoscopic regime 
I > t - t > £ is 



hix^t) = /°(x 1 ,t) + A 1 (x 1 ,t) 

= e- ii ?^)A°(x 1 ,t )+/ t ^e-^^') Jdx^UQrx-r.l) (3.21) 



'to 

(> " Pl °(t'-to), Pl o,to)mr 2 



f°x(ri 



in 



P20 

m 



(t' -t ),p 20 ,to). 
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Note that p 



P'iO- 



The magnitudes of - and - are of course different for different 



-12 



10 



-13 c 



systems. For a dilute gas system, typical values are 10 ~ 10 s and 10 
respectively. The second term of the r.h.s. of ( |3.22| ) is the secular term. Indeed, it 
can be shown that in the spatially homogeneous case it is proportional to t — to [25 . 



Accordingly, we have chosen /i(,to) to be zero following the prescription given in ]2 
The RG equation reads 



0. 



t=t 



(3.22) 



at m or 1 



dx 2 —U(\ ri -r 2 \) ■ — / 1 (r 1 ,p 10 ,t)/ 1 °(r 1 ,p 20 ,t)(3.23) 



' dri 



dPi' 



In ( 3.22j ) we have imposed that t = t although the expression (|3~22|) is valid for t-t > J 
. This manipulation can be justified by the same logic given in the last part in §2 
and will appear also in the case of field theory discussed in the following section: The 
t-derivative is the microscopic derivative and the to-derivative is the macroscopic one. 
Through the RG equation, we can automatically go over to the mesoscopic physics from 
the microscopic physics. Thus the mesoscopic nature of Boltzmann equation is transparent 
in our approach. We have also replaced the argument r 2 of f® with r\ in ( |3.23|) because 
only 1 7*i — 7*2 1 < d contributes to the integration. 

( |3.23|) is the kinetic equation we have been seeking for. In the language of the RG 
method, it is the renormalization group equation describing the slow motion on the invari- 
ant manifold with the coordinate fi(xi,t) [ 23fl . To obtain the usual Boltzmann equation 



which contains the gain minus loss term, we have to manipulate the r.h.s. ignoring the 
spatial dependence. The result is 



0_ 

Ot 



fi(x u t) 



Pi 

m 



d 
dri 

pdp 



2tt 



dp 2 M/°(»"i, Pi> P 2 > t) - /i°(ri, pi, t)fl{r u p 2 , t 

(3. 

where we have introduced the cylindrical coordinate pointing the direction of the relative 
velocity v 12 = (p 2 — Pi)/fn- Comparing (|3.11| ) and ( |3.23|) , we see the change of the 
equation by including the lowest-order contribution of the collision. 



3.2 Role of the initial condition 

Bogoliubov |Q was the first to point out that Boltzmann equation represents the meso- 
scopic physics and derived it from Liouville equation. In his derivation, he assumed that 
the system can be described only in terms of the one-particle distribution function. That 
is, starting from an arbitrary initial condition, the system will rapidly reach the state in 
which f s (s > 2) are functionals of f±. 

f s (xi, ■ ■ -,x s ,t) = f s (xi, ■ ■ ■ ,x s ;fi( ,t)) s>2. (3.25) 
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(f s depends on time only through f x .) In fact, this assumption is a basis of any kinetic 
theory. The special expansion method based on this assumption leads to a special so- 
lution to BBGKY hierarchy. In the RG method, we started with a naive perturbative 
expansion without any knowledge about kinetic theory. Physics enters when we choose a 
special initial condition ( |3.18| ) and with this choice we can construct an invariant manifold 
spanned by the one-particle distribution function in the infinite-dimensional functional 
space, which was originally envisaged by Bogoliubov|Q; see also [03 . 
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4 Derivation of kinetic equation in quantum field the- 
ory 

In this section we apply the RG method to derive Boltzmann equation for quantum field 
theory. In fact, the application of the RG method to quantum field theory has been 



extensively done by Boyanovsky, de Vega and Wang |24| in the prescription given by 
Illinois group ||16||. In the present section, using the self- interacting scalar field theory, one 
of the model Lagrangians which Boyanovsky et al treated, we shall apply the RG method 
as formulated in 



19, 20 



and elucidate how the choice of the initial values is essential to 
derive the time-irreversible equation and that the scale of the time in the RG equation 
is different from the microscopic time in the original equation. Although the formulation 



and the points to be emphasized are shifted from what is given in 24"], we shall closely 



follow the calculation by Boyanovsky et alp4 
Consider the following Hamiltonian: 



H 



H'o 



H'o + H Lv 
1 

2 



Hi 



(fit 



j3 X 

ax— i 
4! 



(4.1) 



(V0) 2 + m^ 



At finite temperature, any particle feels the mean field produced by the thermal excited 
particles which are present even in the classical regime; hence a change of the particle 
picture from the bare one to what we call the quasiparticle. Fortunately, in the case of 
the scalar theory the medium (thermal) effect only affects the mass of the particle. So it 
is convenient to rewrite the Hamiltonian so that the unperturbed one is written in terms 
of the quasiparticle operators as follows; 



H 

Hint 



Ho + H int , 

X - J d 3 x[n 2 + (V0) 2 + m 2 eff( p 2 ] = J d 3 ku k a^k)a(k), 



(4.2) 



+ -5m 2 (j) 2 } 



{x,t) 



7r(x, t) 



k 



d 3 k 

d 3 k 
(2vr) 3 /2 7r 



{k,t)e 



kx 



(k,t)e 



ikx 



7T 



k 2 + ml 



Cff5 



m 2 eif + 5m 1 



(M) 



:m) 



m r 




[k,t) + a\-k,t)}, 



a(k,t) +a T (-fc,t)], 



Here we have shifted the mass of the unperturbed Hamiltonian to the quasiparticle mass 
m e ff which is to be determined self-consistently; in the equilibrium, the result is nothing 
but the effective mass to be given in the hard thermal loop resummation scheme |31 , |24| ; 



d 3 k l + 2n fc 
(2vr) 3 2u k 



k 



k 2 + m 2 



eff' 



(4.3) 
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where nu is the Bose-Einstein distribution function; accordingly, 5m 2 = — |(0 2 (x)). 

For a non-equilibrium situation, the definition of the quasiparticles depends on time. 
We suppose that the distribution function n k (t) = n(t;u> k (t)) of the quasi particles 
varies slowly. Here we have assumed that the system is spatially homogeneous, so the 
evolution equation for the distribution function to be obtained may be rather called the 
rate equation. 

To describe the slow variation, we anticipate that n k (t) is a macroscopic distribution 
function varying slowly like the distribution function P(u, t) governed by Fokker-Planck 
equation discussed in §2; the corresponding counter part of the distribution function 
P(u, t; to) which includes rapid time-variations will be introduced shortly. 

To solve the Heisenberg equation for a(k,t), we impose the following initial condition 
at t = to where to is arbitrary; 

(a\k, t )a(k', to)) = n k (t )5(k - k'), (4.4) 

which implies that the quasi-particles has the mass dependent on the initial time to as 
follows, 

We define the microscopic distribution function around t ~ to by 

n k {t- 1 , [n(to)))6(p) = (a^{k, t)a{k, t)), (4.6) 
with the initial condition 

n k (t ; t , [n(t )]) = n k (t ). (4.7) 

nu(t;to, [n(to)]) is a counterpart of the microscopic distribution function P(x,v,t;t ) in 
contrast to n k (to). In the following, we shall suppress the third argument in n k and write 
as n k (t; to). The equation of motion for the distribution function n k (t; to) is obtained by 
simply differentiating with respect to time and using the Heisenberg equation of motion. 
The expectation value can be evaluated perturbatively by the closed time path formalism 



33], [34]]. To two loop order, the result is 

T n k {t;t ) = - / duR(cu,k; [n(t )\) LV fcA ° ;J , 4i 
at * 3 J-od 7i(u — u k ) 

R(u,k;[n(to)]) = ^- [l[ 7- (27r) 3 ^ 3 (fc - Ql - q 2 - q 3 ) 

2uJ k J i ( 2?r ) 2uJ Qi 



x ^5(u + ujq ± + ujq 2 + ujq^N^to) + 35(uj + uj Qi + ujq 2 - uJq 3 )N 2 (t ) 
+ 35(w - LOq x - uq 2 + uq 3 )N 3 (t Q ) + 5(u - u qi - Uq 2 - Wq s )JV4(t )|, 
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iVi(t ) = (1 + n k (t ))(l + ng i (t ))(l + nq 2 (t ))(l + nq a (t )) - n^(t )nq^t )nq 2 (l^)nq a (to), 

N 2 (t ) = (1 + n fe (t ))(l + n 9l (t ))(l + nq 2 (t ))nq 3 (t ) - n k (to)n qi (t )nq 2 (to)(l + nq 3 (t )), 

N 3 (t ) = (1 + n k (t ))n qi (t )nq 2 (t )(l + nq 3 (t )) - n k (t )(l + n qi (t ))(l + nq 2 (t ))nq 3 (t ), 

N 4 (t ) = (1 + n k (t ))n qi (to)nq 2 (to)nq 3 {to) - n k (t )(l + n 9l (* ))(l + ng a (*o))(l + nq 3 {t )). 

Notice that the distribution functions in r.h.s are all the initial ones. Integration over t 
can be easily performed 

n k (t) = n kM + * f *; W , )])IZi*^^M. ( 4 . 9) 

The second term of r.h.s. can include a secular term. Indeed, owing to the limiting 
behavior familiar in deriving the Fermi's golden rule 

l-cos[(a;-fc)(t-to)] t-to-yo u . u , , , \ (A m , 

— — > (t-to)d(u-k), 4.10 

7r(o; — ky 

^ becomes for t — t — >• oo 

A 2 

n k (t; t ) ~ n k (t ) + — (t — t )R(u k , k; [n(t )]) + non-secular terms. (4-H) 

o 

Notice that the r.h.s. is composed of the macroscopic distribution function. The important 
point is that although we have taken the formal limit t — t — > oo, ( (4.11|) is still a local 
expression valid near t . By "near to" we mean that t — t is small in comparison with 
the macroscopic time scale which is in the same order as the relaxation time. Thus we 
are naturally led to suppose that the time variable to is a macroscopic one. The limit 
t — to — > °o merely means that t — to is large in comparison with the microscopic time 
scale which is in the same order as the inverse of the (thermal) mass of the particle. As 
was emphasized by Boyanovsky et al |24], the separation of these scales is justified in the 
weak coupling regime, and the existence of such an intermediate scale is essential for the 
validity of the kinetic equation like Boltzmann equation. A nice point in the RG method 
lies in the fact that we naturally get to change our standpoint to observe the system from 
the mesoscopic point of view, as noticed above. 

Having made this change of the time scale, the RG equation 



t=t 



(4.12) 



leads to the equation for the initial value nu{t) as the macroscopic variable with the 
macroscopic time t as follows 



n k (t) = T fl H , fc;W( )]) = T ^/np- F |_(2^ (fc _, 1 _, 2 _ 93 ) 

x |<J(w + + Uq 2 + LUq^N^t) + 35(u + u q% + u)q 2 - Uq 3 )N 2 (t) 

+ 3S(ou -ujq i -0Jq 2 +ujq 3 )Nz{t)+5{uj-ujq i -ujq 2 -ujq 3 )N±{t)\, (4.13) 
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which is nothing other than Boltzmann equation governing the slow variation of the 
distribution function. Thus the integration constant nu(t) as given by the initial value 
has become a dynamical variable whose dynamics ( |4.13| ) represents the reduced dynamics 



of the system, in accordance with the general property of the RG method elucidated by 
Ei et al[pCfl. We should remark that Eq. (|4.13|) is supplemented with the time-dependent 
mean-field equation governing the time-evolution of m e g(t), which is now given by 

mlS) = < + \S d ^ l \TS^ W *(<W* 2 + "&(')• ( 4 - 14 ) 

This condition leads to the r.h.s. of (4.8) which only contains the potential secular term 
and is equal to the time- dependent version of (4.3). Thus the quasiparticle defined by the 
hard thermal loop resummation scheme is the relevant slow dynamical variable and we 
must choose the initial condition to be the distribution function of it to obtain the rate 
equation. The importance of the choice of the initial condition is again recognized. 
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5 Fluid dynamical limit of Boltzmann equation 



In this section, we apply the RG method formulated above to obtain the fluid dynamical 
limit of Boltzmann equation |3]]; in other words, we shall derive Euler and Navier-Stokes 
equations successively. 

5.1 Basics of Boltzmann equation 



Boltzmann equation is an evolution equation of the one-particle distribution func- 
tion f(r,v,t) in the phase space, and reads 

Here the left-hand side describes the change due to the canonical equation of motion while 
the right-hand side the change due to collisions; 

/[/] = J dv 1 J dv' J dv[w(v Vx\v' Vi) 

x lf(r, v', t)f(r, v[, t) - f(r, v, t)/(r, v u t)\, (5.2) 



which is called the collision integral. The transition probability w(vvi\v' i/J has the 
following symmetry due to the time-reversal invariance of the microscopic equation of 
motion; 

w(vvi\v' v[) = w(v' v'^Wi). (5.3) 
Furthermore, the invariance under the particle-interchange implies the following equality; 

w(v Vi\v' v[) = w(vi v\v[ v') = w{v\ v'\v\ v), (5-4) 



where (|5.3| ) has been used in the last equality. 

Owing to the conservation of the particle number, the total momentum and the total 
kinetic energy during the collisions, the following equalities hold, 

J dvl[f] =0, J dvvl[f] = 0, J dvv 2 I[f] = 0. (5.5) 

We say that the function (f(v) is a collision invariant if it satisfies the following 
equation; 



dv(p(v)I[f) = 0. (5.6) 

For a collision invariant <f(v), we can define the density and the current j as follows; 

n v = J dvip(v)f(r, v, t), j v = J dvvip(v)f(r, v, t), (5.7) 
which satisfy the continuity or balance equation; 

dtn v + V-j v = 0. (5.8) 

19 



Thus we have formally the following fluid-dynamical equations as the balance equations 
for the conservation of the particle number, total momentum and kinetic energy; 

d t p + V-(pu) = 0, (5.9) 
d t (pUi) + djipujUi) + djPji = 0, (5.10) 



,2 



ft(py + e) + di 



,2 



\p— + em. + Qi 



+di(P ijUj ) = 0, (5.11) 



respectively. Here, p, u, e, P^, Qi are the particle density, the stream velocity, the energy 
density, the pressure tensor and the heat flux, respectively and defined as follows; 



p(r 


t) = m 


p(r,t)u(r 


t) = m 


e{r 


w 


Pij(r 




Qi(r 





vf(r,v,t)=mn(r,t), (5-12) 

J dvvf(r,v,t), (5.13) 

dvj\v-u\ 2 f(r,v,t), (5.14) 

dvm(vi - Ui)(vj - Uj)f(r,v,t), (5.15) 

dv—\v-u\ 2 (vi-Ui)f(r,v,t). (5.16) 

These equations are formal ones in the sense that the distribution function / is not yet 
solved: The solution obtained from Boltzmann equation will give the explicit forms of the 
internal energy, the transport coefficients and so on. 
The H function is defined as follows: 

H(r,t)= Jdvf(r,v,t)(\nf(r,v,t)-1). (5.17) 

For equilibrium states, the H function is equal to the entropy S with the sign changed. 
Defining the corresponding current by 

J H (r,t) = J dvvf(r,v,t)(\nf(r,v,t)-l), (5.18) 

one has the balance equation; 

an ~ 

— + V-J H = Jdvl[f]\nf. (5.19) 

This shows that when In / is a collision invariant, the H function is conserved. In fact, one 
can show that this condition is satisfied when f(r,v,t) is a local equilibrium distribution 
function as given by fl5.32j ) below. 

5.2 Derivation of fluid dynamical equations in the RG method 



To make it clear that the following discussion fits to the general formulation given in |2^ 



we discretize the argument v as v — > Vi^: Discriminating the arguments (r,t) and V\ in 
f(r,Vi,t), we indicate Vi as a subscript i for the distribution function; 

f(r,v l ,t) = f l (r,t) = (f(r,t)) l . (5.20) 
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Then Boltzmann equation reads 



?A = iifh _ Vi . ?A (5.21) 

dt [Ji dr' K J 



where 



I[f\i = E ™(«i«>*«0(/fc/j - fifj){r,t). (5.22) 

Now let us consider a situation where the fluid motion is slow with long wave-lengths 
so that 

Vi -y± = 0{e), (5.23) 

where e is a small quantity, |e| < 1. 

To take into account the smallness of e in the following calculations formally, let us 
introduce the scaled coordinate r defined by 

d d 

•■ = «■. * = (5 - 24) 



Then (|5.21|) now reads 



which has a form to which the perturbation theory given in [2(J is naturally applicable. 

In accordance with the general formulation given in p0 |, we first expand the solution 
as follows; 

Let fi(f,t;t )) be a solution around t ~ t given by a perturbation theory with fi(f,t Q ) 
being the initial value at t = to] 

f i {f,t = t ;t )=f i {r,t ). (5.26) 

We expand fi(r,t;t ) as 

fi(f, t- 1 ) = if (f , t; t ) + (r , t, t Q ) + ■ ■ ■ , (5.27) 
and the respective initial condition is set up as follows, 

jf(f,t = t ; to) = ri l \r,t ), (I = 0, 1,2...). (5.28) 
The 0-th order equation reads 

dfT /rri (o), 



at 



W %■ (5-29) 
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Now we are interested in the slow motion which may be achieved asymptotically as t — ► oo. 
Therefore we take the following stationary solution, 



(0) 



at 



o. 



which is a fixed point of the equation satisfying , 

(/[/ (0) D, = o, 



(5.30) 



(5.31) 



i(o) . 



for arbitrary f. Notice that (5.31) shows that the distribution function / is a function 
of collision invariants. Such a distribution function is a local equilibrium distribution 
function or Maxwellian; 



fT\f,t;t )=n(f,t ) 



in 



2nk B T(r,t ) 



3/2 



exp 



m\Vj — u(r, t 



2nk B T(f,t ) 



(5.32) 



Here, the local density n, local temperature T, local flux u are all dependent on the initial 
time t and the space coordinate r but independent of time t. 
The first-order equation reads 



d_ 

dt 



(i) 



—Vs 



dfl 



(0) 



dr 



Here the linear operator A is defined by 



nf ] \f m ]=t§ 



f=f 



(5.33) 



(5.34) 



The operator A satisfies the following; Let f\ = f^ipi. Then 

{Af)i = ff ] E w{ij\kl)ff\A + *j) = ^\Mh 



(5.35) 



Here A is an operation acting on ifi; ipi = Defining the inner product between tp 

and ij) by 



{<P,*I>) = J dvcpip, 
one can show that A is self-adjoint; 

(<P,Ail>) = {Aip^). 

One can further show that the five invariants m, v, ^v 2 span the kernel of A 



J i "1 2 

KerA = {m, v, ^-v 2 }. 

The other eigenvalues are found to be negative because one can show 

(<P,A<p) < 0. 



(5.36) 



(5.37) 



(5.38) 



(5.39) 
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We write the projection operator to the kernel as P and define Q = 1— P. 
Applying the general formulation in [20|, one can readily obtain the first-order solution, 

(5.40) 



d~f i0) d~f i0) 
-{t-t )Pv- ^— + A~ 1 Qv J 



dr '" dr 

The perturbative solution up to the e order is found to be 

~l(o) 

~f(f , t, t ) = / (0) (r, t, t ) + e[-(t - t )Pv ■ + A^Qv ■ -|— 



Notice the appearance of a secular term. 

If one stops the approximation and apply the RG equation, df/dt \ 



t =t 



(5.41) 



0, one has 



a j(o) oi(o) 
+ ePv ■ 



Ot 



dr 



0. 



(5.42) 



This is a master equation from which equations governing the time evolution of n(r,t), 

u(r, t) and T(r, t) in / , i.e., fluid dynamical equations: In fact, taking an inner product 
between m, mv and mv 2 /2 with this equation, one has fluid dynamical equations as given 
in (|5.9| )- (|5.11|) . The only difference from ( |5.9|) - (|5.11|) lies in the fact that / is now explicitly 
solved and the energy density e, the pressure tensor and the heat flux Qi are given as 
follows; 

e(r,t) = Jdv^\v-u\ 2 f°\r,v,t) = h B T(r,t), (5.43) 
Pij(r,t) = J dvm( Vi - Ui){vj - Uj)f {0) {r,v,t) = nk B T(r,t)5 ij = P(r, t)^-,(5.44) 
Qi(r,t) = J dv r ^\v-u\ 2 {v l -u l )f^\r,v,t) = Q. (5.45) 

We have defined the pressure P using the equation of state for the ideal gas in the second 
line. There is no heat flux because the distribution function f(°> in the formulae is the one 
for the local equilibrium. Inserting these formulae into (|5.9|) - (|5.11|) , we end up with a fluid 
dynamical equation without dissipation, i.e., the Euler equation in this approximation; 



d(pui) d d 



9,2 \ 9 



u 



(p— + e + P) Ul 



0. 
0. 

0. 



(5.46) 
(5.47) 

(5.48) 



We notice that these equations have been obtained from the RG equation (|5.42|) . It 
should be emphasized however that the distribution function obtained in the present 
approximation takes the form 



f(r,v,t) = f^(r,v,t) + A- 1 Qv 
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df^(r,v,t) 
dr 



(5.49) 



which incorporates as a perturbation a distortion from the local equilibrium distribution 
and gives rise to dissipations. 

One can proceed to the second order approximation straightforwardly and obtain fluid 
dynamical equation with dissipations as the RG equation. The perturbation equation in 
the second order reads 

d-A)pA=- Vl .^. (5.50) 



dt J i dr 

Here, we must make an important notice: We have actually used the linearized Boltzmann 



equation[]35| neglecting the second-order term of / in the collision integral: It is known 
that the neglected term produces the so called Burnett terms which are absent in the 
usual Navier-Stokes equation [RO]. 



Inserting the first order solution (|5.4C ) into |5.50|) , we have for the second order solution 



/ (2) (f,t;t ) = -(t - toKA^Qv ■ VPv- V + Pv- VA^Qv- V}/ (0) (r,t 



-{(A-'Qv ■ V) 2 - \{t - t )\Pv • V) 2 }/^(r, t ). (5.51) 

Adding up the all the solutions up to the second order and applying the RG equation 
df(r,v,t;t )/dt \ to=t = 0, we have 

^ + e{PvVf^+A^QvV^} 

+ e 2 {A~ 1 Qv ■ VPv • V + • VA^Qv ■ V}/ (0) = 0. (5.52) 

Notice that the time derivative of the distribution function in the first order terms in e 
must be retained in this order, because the time derivative is at most of order e. Applying 
the projection operators P and Q from the left, we have 



dt 



+ ePv ■ V/ {0) + e 2 Pv ■ VA^Qv ■ V/ (0) = 0, (5.53) 



eA^Qv V{^— + ePv- V/ (0) } = 0, (5.54) 

respectively. Clearly, the second equation follows from the first one in this order. The 
third term in the first equation represents dissipations. Taking inner products between the 
first equation and the particle number, the velocity and the kinetic energy, one will obtain 
a fluid dynamical equation with dissipations included, i.e., Navier-Stokes equation |35[ . 

In summary, we have shown that the fluid dynamical limit of the Boltzmann equation 
can be obtained neatly in the RG method as formulated in [P0|| . Such problems to reduce a 
kinetic equation to a further slower dynamics appear quite often, reflecting the hierarchy 
of the space-time of nature. A typical problem [§, [IT], [4(| in this category 



is to derive Smoluchowski equation|4T| from Kramers equation^]. As is expected, it is 



possible to develop a systematic theory for the adiabatic elimination of fast variables in 
Fokker-Planck equations, which is described in the next section. 
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Adiabatic Elimination of Fast Variables in Fokker- 
Planck Equation 



In this section, we shall show that the RG method as formulated in [2(J gives a systematic 
way for the adiabatic elimination of fast variables appearing in Fokker-Planck equations 
f36| , ||, |37|, |38|, [39| . A typical problem in this category is to reduce the so called Kramers 
equation[^2| to Smoluchowski equation JO]]. Although Brinkman|p6| was the first for a 
serious attempt, reliable derivations was given rather lately || |37|, |38], ^9|; see || (|(J for 
a review. Recently, it was shown that the RG method as formulated in [16| can be used 
to eliminate a fast variable in the Fokker-Planck equation[H3|. We shall show that the 



RG method gives the results which have a clear correspondence with those given in other 
methodsH, H(J; thereby, we would say, our method gives a foundation for these methods. 



6.1 A generic FP equation with a fast variable 

Using a generic example, we shall present the method for adiabatic elimination of fast 
variables based on the RG equation. 

The example is the following 2-dimensional Langevin equation with 7 being a large 
number, 

x = h x (x,y) + g x (x,y)T x (t), 

y = jhy(x,y) + f(x,y) + y/ig y {x,y)T y (t), (6.1) 

where Ti(t) (i = x,y) are independent Gaussian noise satisfying 

(T i (t)r j (t')) = 28 ij 8(t-t'). (6.2) 

This example is treated by Gardiner ||. We follow Gardiner for the notations. We suppose 
that owing to the large friction 7, the variable y is a fast variable. Our task is to eliminate 
the fast variable adiabatically and obtain the reduced dynamics for the slow variable 
x only. Before doing this task, we first eliminate the most rapid variables, i.e., r\(i) 
(i = x, y) to obtain the dynamics written only in terms of x and y. This is tantamount 
to transforming the Langevin equation to Fokker-Planck equation as was done in §2; it 
implies that what we are trying to is to further reduce a kinetic equation to a slower 
dynamics. The corresponding Fokker-Planck equation for the probability W(x, y, t) reads 

dW 

\L 1 + 1 L ]W } (6.3) 



where 



dt 



In = i^P, • ^fU.yn • |U,- (6-4) 
d -D x + -^- 2 D yy ), (6.5) 



with 



' dy dy 2 



Di = hi(x,y) + gi{x,y)—gi(x,y), D u = g 2 (x,y), (i = x,y). (6.6) 
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To apply the perturbation theory, we introduce a scaled time r by 

er = t, (6.7) 

with e = I/7 < 1 and write r as t for the moment. Our equation now reads 

(d t - L )W = eL.W, (6.8) 

We denote the eigenvalues and the eigenfunctions of the unperturbed operator L (y, x) 
as follows; 

L (p n (y; x) = -X n (x)(p n (y; x). (6.9) 

Notice that x in L (y;x) plays a role only as a parameter. The corresponding adjoint 
equation reads 

LWniv, x) = -\ n {x)<fi{y; x). (6.10) 

We assume that X n (x) > and that the 0-eigenvalue A = is non-degenerate: The 
eigenfunction (p (y;x) is proportional to the stationary distribution function, so we can 
write as <po(y;x) = W st (y;x). It can be shown that the eigenfunctions ip n (y;x) and 
(Pn(y; x) with n > 1 are written in terms of <po(y] x) as 



¥n(y;x) = \/ip (y;x)^ n (y;x), <pl(y;x) = tp n (y;x)/^(p (y;x), (6.11) 
where ip n (y; x) is an eigenfunction of the Hermite operator 



H = , = L y/(p (y;x), (6.12) 

vVo(y;z) 

with the eigenvalue X n (x); 

H ip n {y; x) = -X n ip(y; x). (6.13) 
Owing to the Hermiteness of H , the orthonormality and completeness hold; 

J ' ipn(y, x)ip m (y; x)dy = J <p\(y; x)ip m (y; x)dy = 5 nm , (6.14) 

00 00 

J2^n(y;x)ijj n (y';x) = J2 ( fi(y^ x ) l fn(y';x) = 5(y -y'). (6.15) 

n=0 n=0 

For later use, we define the projection operator P to the kernel of L by 

lP V ] { y; X ) = My-,z)J^M*W- (6.16) 

The projection operator to the orthogonal space to the kernel is denoted as Q, 

00 „ 

q = 1 - p = £ My, x) / dy'My'; x). (6.17) 

n=l J 
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We remark that 

vUy;x) = l. (6.18) 

In the following, we assume that 

PUP = 0, (6.19) 

which can be always satisfied by redefinition of L and L\. In terms of the notions 
in quantum field theory, L\ is normal- ordered with respect to the vacuum |0) = tpo] 
(0|Li|0)=0. 

Let us suppose that a solution is given in the perturbation series; W(t, x, y) = Wo(t, x, y)+ 
eW\{t, x, y) + e 2 W2(t, x,y) + . . .. Then following the general scheme of the RG method as 
given in Appendix A of |2"0fl, we first construct the perturbed solution 



W(t, x, y; t ) = W (t, x, y- t Q ) + eW^t, x, y) + e 2 W 2 (t, x, y; t ) + . . . . (6.20) 
with the initial conditions given at arbitrary time t = to; 

W i (t ,x,y',t ) = W i (to,x,y). (6.21) 
Notice that this implies that 

W(t,x,y) = W(t,x,y;t). (6.22) 
Now the equations for W (n — 0, 1, 2, . . .) read 

(d t -L )Wo = 0, 

(d t -L )W n = InWn-!, (n = l,2,...). (6.23) 

The solution to the lowest order equation around the initial time t — to now reads 

oo 

W (t, x, y; t ) = £ C n (x; t )ip n {y; x)e" A "*. (6.24) 

n=0 

Notice that the coefficient function C n (x;to) may depend on the initial time to as we U as 
x. This solution implies that the initial value is taken as 

oo 

W(t ,x,y) = W (to,x,y;to) = £ C n (x; t )<p n (y; x)e~ x ^ . (6.25) 

n=0 

If we are only concerned with the asymptotic (long time) behavior as t — > +oo, however, 
we may keep only the stationary solution in the sum, namely, 

W (t, x, y; t ) = C (x; t )(po(y; x) = W {t , x, y). (6.26) 

Notice the last equality, which gives the functional form for the true solution in the 0-th 
order, holds because of the stationarity of the solution. We remark also that 

PWoit, x, y; to) = W (t, x, y; t ). (6.27) 
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Let us solve Eq.( |6.23| ) order by order. The solution to Eq.( |6.23| ) is formally given by 

W n (t,x,y;t ) = —^L 1 W n - l = —^{P + Q)L 1 W n ^ (6.28) 
o t - L d t - L 

= X n (t,x,y;t ) + Y n (t,x,y;t ), (6.29) 



where 

X n = PW n {t,x,y;t ) = 1 PLiWn-i, 

o t - L 

Y n = QW n (t,x,y;t ) = —^ r QL 1 W n _ 1 , (6.30) 

o t - L 

with n — 1,2, ... To obtain the special solutions to these inhomogeneous equations, one 
must impose an initial condition. We impose the initial condition at t — to so that 
possible P-space components in the perturbed solutions is "renormalized away" to the 
unperturbed solution|TJ5|, |2"0" |; namely, we shall construct special solutions which satisfy 

PW n (to, x, y- 1 ) = X n (t , x, y- t Q ) = 0. (6.31) 

We shall see that this condition is satisfied owing to the appearance of secular terms 
proportional to (t — to) n multiplied to the P operator. 

The special solutions to the perturbative equations ( |6.23|) with the initial condition 
( |6.31| ) may be obtained in a mechanical way by the method given in Appendix A of |20|| : 
For any constant vector U the following formulae hold [^0|] ; 

' -QU = -^-QU, (6.32) 



1 -{t-t ) n PU = 1 (t - t ) n+1 PU, (6.33) 



(6.34) 



for n = 0, 1, 2..., especially for n — 0, 

1 



d t - L 



PU =(t- t )PU, (6.35) 



and 



-{t-t^QLxU = (t-t )^ r QL 1 U - — — QL X U. (6.36) 



Here, notice that [I/o,Q] = 0, so Q/L = L^Q = QL^ l Q. We remark that these special 
solutions are all compatible with the initial condition (|6.31|) . 

In the first order, one immediately sees that X\ identically vanishes because 

PLiWo = PhPWo = 0, (6.37) 
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on account of ( |6.19| ) and ( |6.27| ). Thus we have for the first order solution 

1 



W l (t,x,y;t )=Y 1 



-Li 



QLtWo. 



(6.38) 



A couple of remarks are in order: (i) One may replace the operator Q with P + Q = 1 in 
this expression owing to Eq. fl6.37D . Although such simplification in the expressions may 
be done also for higher-order solutions, we shall retain the Q operator for definiteness, 
since the expression (— I/ ) _1 I/i itself is ill-defined because of the zero eigenvalue of L . 
(ii) The solution ( |6.38| ) satisfies the initial condition ( |6.31| ) as 



Wi(t ,x,y,t ) = W!(t ,x,y) 



QLxW Q (t Q ,x,y). 



(6.39) 



The second order terms are evaluated as follows 

X 2 ■ 



O t — Lq — Lq 

(t-t )PL 1 -^L 1 W , 

— Lq 



(6.40) 



where (|6.35| ) has been used. Similarly, using ( |6.32j ) we have 

Y 2 - 



O t — Lq —Lq 



Q_ 

-Lr 



LtfWQ. 



(6.41) 



Adding these terms, we have Wz(t, x, y, to). We remark that 

W 2 (t Q ,x,y,t ) = W 2 {t ,x,y) = j-QLx J-QLiWo(t , x, y), 

Lq Lq 



(6.42) 



which satisfies the initial condition (|6.31| ). 

Let us proceed to the third order equation. On account of ( |6.19| ), PL\X 2 = 0. So, 



1 



d t - L 



PL!Y 2 , 



Q 



(t-t )PL 1 (^ r L 1 ) 2 WQ J 

—Lq 



(6.43) 



where ( |6.35| ) has been used in the last equality. Similarly using ( |6.32| ) and ( |6.36| ), we have 
Y 3 = [(t-to)-S-^i^i-S-A 

— Lq —Lq 

n . . n . n . _ 

(6.44) 



Q UplQu + ^l^Wq. 



l-Lo) 



~Lq 



-Lq 
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Adding the above terms, we have for the third order solution 

W 3 (t,x,y;t ) = [(t-t ){PL 1 (-2-L 1 ) 2 + -^L 1 PL 1 -2-L 1 } 

—Lq — Lq — L 

Q f n f Q f , , Q f \3 



LlPLl - Ll + ^ Li y ]Woj (6.45) 
(-L r -Lq -L 

which satisfies the initial condition ( |6.31| ). 

Then utilizing the formulae (|6.33| ), X 4 is evaluated to be 

x 4 = [ht-ufPL^L^L^U + it-t^i-pU-S- UpU-^-U 

£ —Lq —Lq {~Lo) ~Lq 

+PL 1 (-2-L 1 ) 3 }]iy . (6.46) 
—Lq 

Similarly with Y 3 , we obtain 

n = [(t-t ){^-^i(^i) 2 + (^i) 2 p^i^i} 

— L — Lq —Lq —Lq 

Q r of , Q ? ^2 Q r Q f , 2nf Q 



-LiPL^^LtY - -^(-^-L i yPL 1 -^L 1 

y — Lo) 2 —Lq —Lq —Lq —Lq 

Q t Q f of <2 f , / <5 f N4n 



L 1 —^-L 1 PL 1 ^ r L 1 + (-^Li) 4 ]^- (6.47) 



-L (-L ) 2 -L -P 

Adding the two terms, we have 

W 4 (t,x,y;t ) = [ht-tofPLx-H-LtP^-H-U 

6 —Lq —Lq 



+(t - tQ){-PL 1T - 2— + PLiC-S-Li) 

+ ^-L 1 PL 1 (^-L 1 ) 2 + (-2-^)^-2-^} 

Q f n f / Q f A 2 Q / Q M2? nf Q 



-LiPLi(-^-Li) -t— (— -t— Li) LiPLi— Li 

v — -^o) 2 — -^o — -^o — -^o — -^o 

<9 f <2 f n f Q f , , Q f v4i 



p - LlPLl - p 1 + ( - Ll f]w Q . (6.48) 

( — -^o) — -^0 — -^0 —Lq 

We remark that PW / 4 (to, x, y; t Q ) = 0. 

If we stop at this order, the perturbed solution is given by the sum 

4 

W{t, x, y; t ) * E ^{t, x, y; t ). (6.49) 

8=0 

This solution becomes invalid a.s t — tQ —>■ oo due to the secular terms. 

Now the RG equation can now re-sum the secular terms: The RG equation 



dW 



dtn 



= (6.50) 

t =t 
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gives for the reduced equations in the P- and Q-space 

3 



d t PW (t, x, y) = e'PL.iY: ^(-S-Li)" - e 2 -S— ^-^Pl^L^W^ x, y) {6.51) 

n=l —-^0 —Lq —Lq 



and 



Q t i .2/ Q f \2i q tj/ _ r,3 Q r or Q ? , ,4 <5 f or / Q r \2 



{e-^Li + e^-^L!) 2 }^^ = {e 6 ^ L X P L x -^ L x + e 4 -^L 1 PL 1 (-^L 1 

— Lq — — — -^0 — -^0 — -^0 



Q , ,2 i.i Q 



+e\^ r L l YPL l ^ r L l }W 0} (6.52) 

respectively. Here we have utilized the fact that W (t,x, y;t = t) = W (t,x,y) and 
Eq.( |6.27|) . It is readily verified that Eq. (|6.52|) follows from Eq. (|6.51|) . Eq. ( |6.51|) is exactly 



the same as the one given by Gardiner || where the Laplace transformation was utilized. 
The fourth-order term in (|6.51|) may be rearranged as a wave-function renormalization in 
this order of approximation; 

(l + e 2 P^-^^)9 i H/o(t,x,y) = e 2 PL 1 [^ e "- 1 (^-L i r]H/ (t,x, 1/ ). (6.53) 

{~^0) n=l —^0 

In fact, multiplying Eq.( |6.52| ) by —ePLiL^ 1 and summing the result with Eq.( |6.51| ), one 
arrives at Eq. (|6.53 ). 



Eq.( |6.51| ) or ( |6.53| ) is actually the equation for the amplitude C (x,t): Performing the 
projection P to the both sides, which is nothing but performing the integration / dy of 
the both sides, one has after recovering the original time t — > t/e with e = I/7, 

(l + 1 ^^ N \x))d t C (x,t) = j2l Hn - 1] £n(x)C (x,t), (6.54) 



n=2 

with 

00 1 

3 N \x) = £4 n -_4 )0 , (6.55) 

n=l ^n\ X ) 

&*( x ) = J2 £o,n^-r-r£n,o, (6.56) 

00 . 1 1 

n,m / —r£>m,0} (6.57) 

n ^ =1 Xn{x) Xm(x) 

00 1 1 1 

n,m \ / \ £"ml \ / \ ^£,0- (6.58) 

„^ =1 ' A n (x) A m (x) A,(x) 
Here L n ^ m (x) are the y-averaged operators with respect to a;, defined by 

C n , m (x) = / dy<fl i (y;x)L x (x,y)ip m (y;x). (6.59) 
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6.2 Derivation of Smoluchowski equation from a Kramers equa- 
tion 



As an example, let us take the following simple case||43||; 

h x (x,y) = y, g x (x,y)=0, h y (x,y) = -y, 
f(x,y) = F(x), g y {x,y) = Jr(x). (6.60) 

This equation is a Langevin equation with a multiplicative noise. In this case, we have 

L = d y y + T(x)d 2 y , L x = ~yd x - f(x)d y . (6.61) 

The eigenvalue problem for L is easily solved; the eigenf unctions belonging to the eigen- 
value X n = n, (n — 0, 1, 2 . . .) may be written 

<p (y;x) = e-y 2 ' 2T ^/Jmx), ip n > x {y;x) = ^H n (y/^{xj)ip (y,x), (6-62) 

where H n (x) is the Hermite polynomial in the n-th order. The conjugate eigenfunctions 
are found to be 

^(y; x) = H n (y/y/2T(x)), rf(y; x) = 1. (6.63) 
Accordingly, the projection operator P is given by a simple integration, 

P = cp (y;x) J dy. (6.64) 

To obtain the reduced dynamics, we only have to calculate the y-averaged operators 
C nj o(x) and Co >n (x). First we notice that PL\P = 0, hence C\ = Cq^x) = 0. The 
non- vanishing term is evaluated to be 



Co,n( x ) = I dyLxy n {y;x), 

dy(yd x + f{x)d y )<p n (y; x), 



-^IVt (6 ' 65) 



So we only have to calculate £ 10 : 



= - I dy(p\(y;x)(yd x + f(x)dy)(p (y;x). (6.66) 
A straightforward calculation gives 

l~2 d dT 
(yd x + f(x)d y )<f Q (y;x) = ip x {y; x) J - {T— - f(x) + — } 

+^ 3 (2/;s)^/P^, (6.67) 
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hence 



4oW=-/|{r|-/W + f}. (6.68) 
Using the above results, one obtains in the second order approximation 

£2 = -^0,1-^1,0) 

dx ( dx ^ ^ ~^ dx ) ^ ^ 



Thus we end up with [|43l 

dC (x,t) 
dt 



r l t 2 c,{x-t), 

= ^fx{ T lx- f ^ + d I) C ^- ^ 

One can proceed to higher orders. The third order correction is found to vanish owing 
to the parity conservation Ci^ = £i >3 = 0. So one must calculate the fourth order to have 
a correction to the second order result. Although it is straightforward to calculate £4 with 
x-dependent T(x), we shall, in the following, take a simple case where T is independent 
of x. To obtain the fourth order correction, one needs to evaluate the following operators: 

It d 

C v = \!<*4-*u<-^)>, (6-72) 
4o = 0- (6-73) 



Thus one has 



dx 2 dx 

which is exactly the same as given in other methods 



A = ^(T|--/) 2 , (6.74) 
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7 Summary and Concluding Remarks 



In this paper, we have applied the so-called renormalization group (RG) method to de- 
rive and reduce kinetic equations; the equations treated include Boltzmann equation, the 
fluid dynamical equation, Fokker-Planck equation for classical dynamics and also the rate 
equation in quantum field theory. In contrast to previous works ][23| , |24jl , our main purpose 
was to elucidate the general structure of the reduction of the dynamical equations in the 
hierarchy of the evolution equations. We have noticed that the significance of the choice of 
the initial value on the attractive manifold which is also an invariant manifold in deriv- 
ing kinetic equations is fully recognized and emphasized by BogoliubovQ], Lebowitzjl^l . 



KubojriJ and Kawasaki ]|14T| , for instance. The notion of coarse-grained time derivative 
was also noticed by Mori|| and others 0, 11|. Our point was that these basic ingredi- 



ents naturally appear in the RG-theoretical derivation of kinetic equations when properly 
formulated so as to respect the role played by the initial condition as formulated in [^, |2(| . 

We have also shown that the further reduction of kinetic equations can be performed 
in a unified manner in the the RG method as formulated in ||20|| . Such problems include 



obtaining the fluid dynamical limit of Boltzmann equation and deriving Smoluchowski 
equation from Kramers equation. 

It is well known that the RG in quantum field theory(QFT) and statistical physics fl5 



21. 22] works well as a powerful tool for obtaining the infrared effective theories with fewer 
degrees of freedom than in the original Lagrangian relevant in the high-energy region. This 
is a kind of the reduction of the dynamics [pt^ . So one could imagine that the RG may 
be applied to derive and reduce kinetic equations in a unified manner. Now we know 
that it is the case. Conversely, does the coarse graining of time as shown in the present 
work appear and play a role in the RG equations in QFT and statistical mechanics? In 
this respect, one may notice the resemblance of the exact RG equation like Wilsonian or 
Wegner-Houghton RG equation with Fokker-Planck equation; the quantum filed theory 
with high-energy degrees of freedom remained corresponds to Langevin equation and the 
exact RG equation to Fokker-Planck. Such an analogy is not fully pursued yet [HI] . 
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Appendix A Another calculational procedure for 
reduction of Langevin to Fokker-Planck equation 



In this Appendix, we shall rederive Fokker-Planck equation for (|2.1|) starting from the 



stochastic Liouville equation in a more elementary way sketched in |23[ for far simpler 
equation; we shall work out a detailed derivation to show how the identification of the 
initial condition is important for obtaining the averaged equation, thereby exhibit clearly 
the similarity of the discussion with the other problems considered in the text. 
We first make a change of independent variables |23| for ( [2.4|) (t, u) — > (r, x) by 



r = t, x = e{u- I dsh(u(s))}, (A.l) 



where e is supposed to be small. Then Eq. ( [2.4|) is converted to 

df/dr = -e[V(gRf) + V • hf], (A.2) 

where V = ^d/dxi. We now try to solve Eq.( |A.2| ) around t ~ Vto by the perturba- 
tion theory. We suppose that the initial distribution f(u,t ) is given at t = to- The 
corresponding solution is written as f(u,t;t ) and is expanded as 

f = fo + eh + e 2 f 2 + ---. 

The initial distribution is also expanded; 

f = k + th + e 2 f 2 + ---. 

The equations in the first few orders read 

dfo/dr = 0, (A.3) 
dfi/dr = -[V(gRf )+V-hf ], (A.4) 
dh/dr = -[V(^A)+V-/iA]. (A.5) 

The first order solution is a stationary one; 

f (x,t;t ) = f {x,t ), (A.6) 

where the initial distribution fo(x, to) is not yet specified. Owing to the absence of explicit 
time-dependence of jo, the first and the second order equations can be readily solved; 

ft 



/i(aj,t;to) = "/ dsV-(gRfo)-(t-t )(V-h)fo, (A.7) 

Jt 

Ux,t;t Q ) = f l d Sl r ds 2 L 1 {s 1 )L 1 {s 2 )fo + lv -hV -hfo(t-to) 2 

J to J to Z 

+terms linear in R. ( A . 8 ) 
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Here Li(s) is defined in Eq.(2.9) in the text. The averaged distribution function P(x, t; t ) 
is now given by 

P(u,t;t ) = (f(u,t;t )), 

= f (u,t ) - e(t - t )(V • hfo(u,t ) 

+e 2 f dsi f ds 2 (L 1 (s 1 )L 1 (s 2 ))fo(u,to) 

+i(t-t ) 2 V hV -hf (u,t ). (A.9) 

For steady noises, the correlation // rfsi ds 2 (L 1 (s 1 )L 1 (s 2 )) can be further reduced as 
was done in §2; 

f d Sl f 1 ds^L^LM) = (t- t )G(t - t ). (A.10) 

Jto Jto 

The RG equation dP / dt \ to= t = gives 

df (u, t)/dt + eV • hf - e 2 G(0)f (u, t )) = 0. (A.ll) 

Recovering the original variables, we finally have 

dfo/dt + V(-hfo) - G(0)/ o = 0. (A.12) 

This is the desired Fokker-Planck equation, provided that the initial distribution fo(u, t ) 
is identified with the averaged distribution function P(u,t). This means that the initial 
distribution function at an arbitrary time t = to before averaging must coincide with 
the averaged distribution to be determined. The initial value may be considered as the 
integral constant in the unperturbed equation, which would move slowly being governed 
by the RG equation. In other words, the averaging is automatically made by the RG 
method. 
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